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Abstract 

In the framework of field approach with the finite time interval the particle oscillations 
in a medium are considered. The absorption in a final states do not lead to dramatic 
suppression of ab transitions. Also we touch on the problem of infrared singularities. The 
approach under study is infrared- free. It is shown that "correction" published recently 
[A.Gal, Phys. Rev. C61 (2000) 028201] is clearly wrong. 
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1 Introduction 



The theory of ab oscillations [1] are based on single-particle model. The interaction of particles 
a and b with the matter is described by potentials Ua,b- Imllb is responsible for loss of 6-particle 
intensity. The wave functions ^^^5 are given by equations of motion. The index of refraction, 
the forward scattering amplitude /(O) and potential are related to each other, so later on the 
standard approach is referred to as potential model. 

In some instances there is a need to consider the ab conversion in the matter followed by 
reaction 

(a — medium) {b — medium) b + c ^ f. (1) 

Here c is the particle of medium (one should sum over all c-particles) ; b + c ^ f represents 
the reaction. The whole process {ab transition, 6-medium interaction) takes place in the same 
layer of matter. 6-particle absorption is essential and Imllb cannot be ignored. An example is 
the nn transitions [2] in the medium followed by annihilation 

{n — medium) {n — medium) f, (2) 

where / are the annihilation products which should be detected (see Fig. la). Standard cal- 
culation (see, for example [3,4] and Section 2) predicts a dramatic suppression of this process 
due to n-medium interaction. Below it is shown that process (2) is qualitatively equivalent to 
free-space process given in Fig. 2: 

n ^ n ^ pe^v. (3) 

Therefore, process (3) should be drastically suppressed by the decay in a final state, what is 
incorrect. We cite a body of other arguments (see Sec. 6) which point to the fact that potential 
approach is reasonable as starting point only. 

In the framework of field approach with finite time interval (FTA) we perform the direct 
calculation of the processes shown in Fig.l. The connection between FTA and ^'-matrix theory 
is studied. The approach reproduces all the results in neutrino oscillations in which lmfb{0) is 
ignored as well as standard calculation of the process (2). However, we disagree with standard 
calculation, because, in our opinion, the absorption is described improperly in this case. 

2 Potential model 

In the standard approach the nn transitions in the medium are described by Schrodinger 
equations 

(idt- Ho)^ri(x) ^e^n(x), 
{idt-Ho-V)^^{x)^e^rr{x), 
Ho = -VV2m + Un, 
V^Un-Un^ ReUn - ir/2 - Un, (4) 

e = (m2 — mi)/2 = l/Tnn. Here m and Un = const are the neutron mass and potential, 
respectively; Un and F ~ 100 MeV are the optical potential and annihilation width of n; Tnn 

is a free-space nn oscillation time; m\^2 ^-re the masses of the stationary states 77-1^2 • e is a 
small parameter. The initial conditions are ^„(0,x) = n(0,x), ^n(O) = 0, where n{x) — 
y~^/2gxp(-ie„t + ipnx), en = P^/2m + Un- 
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For analogy with field approach we introduce the evolution operator U{t) = I + iT{t). For 
V = const, in the lowest order in e we obtain matrix elements Uii{t) =< n(0) | ^'n(^) > and 
Tni{t) =<n{0)\^n{t)>: 

Tn.{t) = {e/V){exp{-tVt)-l), 
Tu{t) = i{e/Vf[l - iVt - e^v{-iVt)]. (5) 

The probabihty to find antineutron Wn{t) is 

Wn{t) =|T^,(t) |'= (e/ \V\f[l -2cos{ReVt)e-^"^ + e"^*]. (6) 

The unitarity condition gives 

l^Uu{t)\^ +W^{t)+Wann{t), (7) 

where Wannif) is the probability to find the annihilation products. In the potential model the 
probability of nn conversion Wpot is defined by equation 

Wpot{t) = 1- I Uu{t) |'= 2ImTii{t). (8) 

For brevity, we put below 

ri>i. (9) 

Then 

Wpotit) « ^eH/r. (10) 

Noting that free-space nn transition probability is 

Wf = eH' (11) 

for suppression factor Rpot wc have Rpot = Wpot/Wf = A/Tt <^ 1. 

We would hke to stress that Eq.(7) is automatically fulfilled only for correct model with 
Hermitian operators. Also recall that instead of realistic operator of n-medium interaction Un 
the parametrization 

Un = ReUn - ir/2 = const (12) 

is used. On the reasons given below we abandon the approximation (12). In this case the 
usual procedure connected with diagonalization of mass matrix cannot be realized. However, 
the result do not depend on the basis. All existing calculations have been done in n — n 
representation. 

3 Field approach 

Let us consider the process (2) in the framework of field approach. The interaction Hamiltonian 
involves two terms: 

Hnnit) = e J d^xi^n'^n + H.C.), 
H{t) = (all n — medium interactions) — Un, (13) 

Hj = Hnn + H . Here \t'„ and \I'^ are the fields of n and n, respectively; = m„ = m. The 
background nuclear matter field Un is included in unperturbed Hamiltonian Hq; quadratic term 
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Hnn included in Hi. The sole physical distinction with model (4) lies in description of n-medium 
interaction H. If it is putted that H{t) = J (f x{—iT /2)'^ n'if n (what is at least unjustified), then 
the potential model results are reproduced (see below) . Field approach allows one to calculate 
directly the off-diagonal terms. 

Let us Un — ReUfi = 0, F = F^, H = Hp, where F^ and Hp arc the width and Hamiltonian 
of free-space decay n — > pe^v. In this case Eqs.(4) and (13) describe the free-space process (3) 
in the framework of single-particle model and field approach, respectively. (As with process 
(2), the neutron decay is excluded. In this connection the free-space process (3) is imaginary 
one, which is inessential.) Then according to Eq.(lO) the process (3) is drastically suppressed: 
Wpot{t) ~ A^eH/Vp. So in the field approach this process (see Fig.2) should be also suppressed 
by the decay in a final state, what seems unrealistic. 

For Fig. la (and, therefore, for Fig.2) it was shown [5] that the process probability is W{t) ~ 
Wf{t). Below this result is derived from consideration of more general process shown in Fig. lb. 
However, first of all we verify FTA, namely, reproduce the potential model results (5). 



4 Calculation 

First of all we consider the nn transitions with n in the final states (n are detected). The 
similar problem takes place in neutrino oscillations. Due to the zero momentum transfer in the 
vertex corresponding to Hnn{t) the process amplitude is singular 

where M is the amplitude of n-medium interaction. (The amplitude of potential model is con- 
sidered in Sec. 5.) For solving the problem the approach with finite time interval was proposed 
[4] . In the lowest order in e we have 

<riO| U{t)-I\On>^iTf,i{t) = {-i) <npO\ dtpH^fiitp) +T''{t - 0) J^' dtpHr,n{tp) \Onp>, 
T%t-tp)^Y.i-^f dti... ~ dtkH(t^)...H(tk), 

k=l •^*/3 -^^/S 

where | Oup > and | Ofip > are the states of the medium containing the neutron and antineutron 
with 4- momenta p = (e„, Pn). Taking into account that Hnn \ On^ >= e | Oup >, we change the 
order of integration [4] and obtain 

Tni{t) = -et - eJ^dtpiTtiit - tp), (14) 
iTtiir) ^<npO\T^ir)\Onp>, 

where t = t — tp. The n-medium interaction is separeted out in block T^{t). 

For verification of FTA we calculate in the framework of potential model: V — const, 
H{t) = V{t) = exp(i//o^)l^exp(-i//o^) = V. We have 

tTJiiT) = U^ir) - 1 = eM-^VT) - 1. (15) 

Substituting this expression in Eq.(14) one obtains Eq.(6). This means that FTA should re- 
produce all the results in neutrino oscillations. (Commonly, in neutrino oscillations lmfb{0) 
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is ignored.) One further important test of FTA (calculation of Tii{t)) have been given in [5]. 
Therefore, the FTA was verified by the example of exactly solvable potential model. It is 
involved in <npO \ T^\t) \ Onp> as a special case. 

Let us consider the process (2) wherein annihilation products are detected, namely, the nn 
transitions in the nuclear matter (Fig. la). Here, F ~ 100 MeV and n inevitably annihilates. 
We consider the more general problem. We calculate Fig. lb on the interval (t/2, —t/2). As a 
result it will be shown that: (a) When g — > (g is 4-momenta of particle escaped in the nn 
transition vertex), the result converts to one corresponding to Fig. la, which is interesting for 
us. (b) When g 7^ and i — > 00, the result coincides with ^'-matrix one. Such scheme allows 
to verify and study the FTA. (c) The functional structure of Wpot{t) is wrong. 

Let us consider the imaginary process 

n^ri + $. (16) 

For decay to be permissible in vacuum put = m— 2m$. The corresponding process in nuclear 
matter is shown in Fig. lb. This is a nearest analogy to the process under study. Instead of 
Eqs.(13) we have Hi = H'^^ + H, 

Hnnit) = ^' j d?x{^>^^*^>r, + H.c.), (17) 

where if^^ is the Hamiltonian corresponding to decay (16). (For Fig. la H'^^ — > Hnn-) In a 
manner like the calculation of T'm(^) we have 

rtk 



n{t) = I T^{t) I \ dt^KMe-^y \Onp>, (18) 

00 „t/2 ft^-i ft/2 

T^it) = E(-0' / dh... / dtkH{ti)...H{tk) = Texp(-i / dtiH{h)) - 1. (19) 

^ J -t/2 J -t/2 J -t/2 

Here < / | represents the annihilation products with (n) mesons. Multipliere exp(— a 1 1/3 |), 
q; > is introduced for realization of adiabatic hypothesis. One obtains 

T,{t) = le' V < / I r"(i) I Onp_, > NF, (20) 
/\q — la 

p ^ g-a|tfc| _ g-at/2g-iAg(tfc+t/2)^ ^21) 

^9 = 9o ^ 2m$ + (p„ — q)^/2mfj — p^/2m. Here | Onp_g > is the state of the medium containing 
the n with 4-momenta p — g. The normalizing factors of the wave functions of n and annihilation 
mesons are included in < / | T"(t) | Onp_q > and the other those in the multiplier N. 

Nonsingular diagram. - Since Ag 7^ the limit t — > 00 can be considered. Then F — 
exp(— a \ tk\) and is the usual S'-matrix element. FTA reproduces the S'-matrix result, as we 
set out to prove. 

It is easy to estimate the widths corresponding to Fig. lb and free-space decay (16): 

^free ~ e'^m$/(27r), 

where we have put m^/m <^ 1. The ^-dependence is determined by exponential decay law 

Wl f = 1 — p-^bjreet ^ F, , f 

''bjree — J- c ± bjree^- 
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These formulas will be needed below. 

Singular diagram. - Let us return to T(,(t) and consider the formal limit g — > 0. For nn 
transition Hamiltonian the adiabatic hypothesis should not be used: a — (see below). Now 
factor F/Aq is 

[1 - (1 - iAq(t/2 + tk))]/Aq = i / ' dt^. 

J~t/2 

Changing the integration order [4] and going to interval (t, 0) one obtains 

nit) = -e'N f dtp<f\ T'^it - tp) I Onp > . (22) 

J 

At a point ? = 0, i?^^(t) = i^nn(^) and = 1. We have 

n{t)^Ta{t) = -e l^driT^r), 

iT%{T)^<f\T%T)\Qn,> . (23) 

TJ^{t) is an exact annihilation amplitude. Ta{t) coincides with the second term of Eq.(14) with 
the replacement <i |=<fipO |— 1. This can be considered as a test for Tait). (In the Ref.[5] 
only Fig. la was calculated. Also it was shown that there is a double counting in the potential 
model.) 

5 Infrared singularities and S-matrix problem formula- 
tion 

The amplitude of potential model obtained by means of (S'-matrix approach is not singular. 
Prom microscopic theory standpoint the reason is as follows. Due to zero momentum transfer 
in the e- vertex this amplitude contains singular propagator of n. However, it also contains 
block which is a sum of zero angle rescattering diagrams of n. As a result the self-energy 
part E = V" appears (see Eq.(22) of Ref.[5]). We are interesting in off-diagonal matrix elements 
which do not contain above mentioned sum {TJ^ instead of T^^) and hence diverges, because one 
singular propagator after e-vertex appears in any case. (Note, that formal sum of series over e 
gives meaningless sclf-cncrgy part E ~ 1/0.) 

FTA is infrared-free. It naturally connected with conditions of experiment. Really, measure- 
ment of any process corresponds to some interval r. So it is necessary to calculate Ufiir). The 
replacement U (r) — > S'(oo) is justified when the main contribution gives some region A < r, 
so that 

Ufi{T > A) = Ufi{oo) = Ufi = Sfi = canst. (24) 

The expressions of this type are the basis for all ^S-matrix calculations. In principle, the three 

cases are possible. 

1. There is bound to be asymptotic regime, however it is not achieved automatically. Then 
the adiabatic hypothesis is used and usual scheme realizing in the field theory or non-stationary 
theory of scattering takes place. Fig. lb corresponds to this case. 

2. There is no asymptotic regime. Example is provided by oscillation Hamiltonian Hnn{t). 
S'-matrix approach is inapplicable. We have usual non-stationary problem. Because of this, for 
Fig. la the result has been obtained in the framework of FTA. In the vertex corresponding to 
Hnnit) the adiabatic hypothesis was not used. 



6 



3. For the interaction Hamiltonian Hj{t) the asymptotic is reached automatically without 
resort to adiabatic hypothesis. An example of this type is T^{t): T^{t ^ l/F) = i. 
For non-diagonal matrix elements the asymptotic should be reached as well: 

mr > l/F) = Tf, = const. (25) 

Really, annihilation of n in nuclear matter can be considered as a system decay with the 
characteristic time ~ 10~^^ s. On the other hand the observation time t = Tq ~ 1 yr [6] is far 
more then that in the experiments on particle decays. So in our case expression (24) is more 
obvious then for any free-space reaction or decay. Recall that all calculations in the particle 
physics based on condition (24). (For process (3) it has a form T^- = const, where is matrix 
of /3-decay.) Besides, for it was verified by direct calculation. 
Since t » l/F, by means of (25) one obtains 

Tait) ^ -ietT%. (26) 

TJ,' can be calculated in the framework of S'-matrix theory. However, in our case we know that 
the n- nucleus decay probability is PF" = Y,f^i \ TJ^ P= 1- Finally 

Wann{t) ^ E I -^^^f, eH^W = eH^ = Wf. (27) 

m 

Due to the annihilation channel nn conversion is practically unaffected by the medium. So 
Tnn ~ Tnn, whcrc T„fj is the oscillation time of neutron bound in a nucleus. The t-dcpcndence 
of Wann IS determined by that of more slow subprocess of the nn conversion Wf ~ t^. Formally, 
quadratic dependence follows from expression for Ta{t). 

All the results have been obtained by means of formal expansions. They are valid at 
any finite t. Consequently, the singularity of S'-matrix amplitude is a direct result of 
incorrectness of the formulation of the S-matrix problem. If t ^ oo Eq.(26) diverges just as Mg 
does. So infrared singularities point to the fact that there is no asymptotic regime. Conversely, 
the analysis of condition (24) suggests the method for handling. The lack of asymptotic leads 
to FTA which is infrared-free. (There is no asymptotic regime for free-space X^X'^-oscillations 
as well. In our opinion, it makes sense to look at the calculation of Am — — mg (GIM 
mechanism) from this standpoint.) 

6 Discussion 

From Eq.(lO) it is seen that: 

(1) Wpot ~ l/F, whereas Fig.l gives the inverse dependence | T^^b p~| TJ^ p (see Eqs.(20), 
(26)). Such structure is typical and for f ^ i uniquely determined. 

(2) The t-dependence of the process probability in the medium and vacuum is different: 
Wpot ~ t, Wf ~ t^. It is beyond reason to such fundamental change. In our calculation the 
t-dependence is the same: when q = (Fig. la), Wann ~ t"^ and Wf ~ t^; when q ^ (Fig.lb), 

W},'^ t and Wfree ^ t. 

(3) In the medium the nn transition suppressed by a factor Rp^t — Wpot/Wf — A/Ft. If 
t = Tq = 1.3 yr [6] (Tq is observation time in proton-decay type experiment) and F = 100 MeV 
then Rpot ~ 10~^^ which seems absolutely unreal. The distribution (27) gives i?a ~ 1. For 
related problems typical suppression factor is i? ~ 1. The process (16) suppressed by a factor 



7 



Rb = Th/Tfree ~ l/vT, where the value m$ = F ~ 100 MeV was used. The reahstic example is 
a pion production pn — > ppvr" in vacuum and on neutron bound in nucleus. When is in the 
region of resonance, the pion absorption due to interaction in the final state is very strong. This 
effects on the number of pions emitted from the nucleus, but not on the fact of pion formation 
inside nucleus, which is interesting for us. 

(4) By means of diagram technique it is easy to get the S-matrix amplitude of potential 
model Mpot- From the expression for Mp^t (see Eq.(22) of Ref.[5]) it is seen that there is a 
double counting with respect to H. 

Points (l)-(4) suggest that the functional structure of Wpot is wrong. When annihilation 
decay probability of n-nucleus W"'{t) = Y^f^i \ TJ^ij) 1, the error increases up to 100 
percent resulting in change of functional structure: Wann ~ ^ Wpot ~ t. As W^{t) 0, 
that is ImUn 0, the error also tends to zero and Eqs.(4) give an exact result. 

Solving Eqs. (4) by method of Green functions one can get 

Tii{t)^i^ f dh r dTUf,{T), 
Jo JO 

n^t) = -e fdrU^ir), (28) 
Jo 

where UII{t) is defined by Eq.(15). It is seen that the main contribution gives the region 
T < 2/F. However, in this region it is meaningless to speak about potential Un. 

(The sclf-cncrgy part H = V is due to of zero angle rcscattcring of n. Such scheme is 
artificial, because cr^^ > cr";^ and in the first act of n-medium interaction annihilation takes 
place. Also recall that for sufficiently early times the exponential decay law is violated [7].) 

One might hope that correct result can be provided by appropriate parametrization. How- 
ever, Ufi has been fitted to radically different problem. Really, the coupled Eqs. (4) give rise to 
the following equation: 

{d^ + idt{V + 2Ho) - Hi - HoV + e')*„ = 0. (29) 

is suffice to get Wpot- Meantime, Un is fitted to p-atom and low energy scattering data. 
These problems are described by stationary equations of Schrodinger type 

(VV2m + iT/2)^n = -E^n- (30) 

The distinctions between Eqs. (29) and (30) are obvious. Eq.(29) is even not Schrodinger type. 
It describes n rather then n. For Eq.(29) the papameter Un is uncertain; the physical sense of 
ImUn is not clear. The F-dependence of the results is inverse: Wpot ~ 1/r, whereas Eq.(30) 
gives W ~ 1 — cxp(— Fi). The structure /q drUlKr) (see Eqs. (28)) appears only in the oscillation 
problem. It is necessary to know UII{t) at a scale r ~ 1/F. So it is beyond reason to use Un in 
the Eqs. (29), (4), because it was fitted to absolutely different problem. 

Let us explain this point. Fitting potential Un we fit the matrix elements associated with 
Eq.(30). However, these matrix elements differ radically from those involved in Eqs. (7), (8). For 
example, in the Born approximation the zero angle p- nucleus scattering amplitude is /(O) ~ 
/ d^xUn{:x.). (Compare with Eqs. (28).) Thus from a pure phenomenological standpoint Tjj(t) 
is uncertain as well. It can be calculated only by means of correct model with Hermitian 
operators. 

For n in nuclear matter the f-dependent equation is 

{idt-m + iT/2)^n = ^, (31) 
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^fi.(O) = 1. From the condition of probability conservation we have Wann = 1~ I P= 
1 — exp(— Ft). F is extracted from experiment. When F increases, Wann increases as well; so 
it will also be for oscillation problem. The potential model gives inverse tendency (10). In the 
problem corresponding to Eq.(31) the value of F affects the branching ratio of channels only. 

Taking into account above given analogy with process (3) we conclude: when ImUn 7^ 0, 
Wpot is wrong. The scheme based on Eqs.(12),(7) is improper. The parametrization Un ~ —iT/2 
is very useful for Eqs.(30),(31) i.e., for the problems with prepared n. However, it is incompatible 
with a realistic n-nuclear dynamic characteristic of nn conversion. 

Finally, we point to an important difference between calculations. Since the asymptotic of 
UII{t) is C/j"(t 3> 1/F) — >• 0, the values of Ta and Wpot are defined by behaviour of UII{t) at 
a scale r ~ (10~^^ — 10~^^) s, which is very "undesirably". For off-diagonal matrix elements 
the picture is inverse: TJ^^ij ^ 1/F) — >• const 7^ 0. (Otherwise, the corresponding process 
probability is VFjj(T) =| T'j\{t) p= 0.) As a result the main contribution to Eq.(23) gives the 
region 1/V < t <t (rather then r < 2/F), resulting in Eqs.(26),(27). Considering off-diagonal 
matrix elements we obviate the principal difficulties mentioned above. 

In fact the problem is very original. From a view-point of optical potential it is that we 
deal with the system of time-dependent coupled equations. In consequence of this the results 
are expressed through the t-dependent matrix elements, which should be calculated beyond the 
potential model. Hence we come to the field approach. However, the oscillation Hamiltonian 
corresponds to 2-tail (because of this the single-particle model is used) which leads to infrared 
singularities. To avoid them the problem is considered on the finite time interval. So we come 
to the time-dependent description again. 

7 Comment on Gal's paper 

In the paper [8] our result was improperly "corrected". The probability to find the annihilation 
products is Wann{t?) ^ eH'^W^ (see Eq.(27)), where is the n- nucleus decay probability 
PF" = 1 — exp(— Ft) ^ 1; e^t^ is the free-space nn transition probability. "Correction" is that 
instead of W"" the probability to find n (Eq.(6)) is substituted. As a result Wann{t) ~ ^"^t^Wn. 
(See Eqs.(31)-(33) of Ref. [8].) We read: The probability to find annihilation products=(the nn 
transition probability) x (the probability to find n-nucleus). Obviously, this "correction" is only 
trivial error. The rest of the "results" [8] are known since 1980. For example, Eq.(27) of Ref. [8] 
(the final result) coincides with Eq.(2) of Ref. [5] (the beginning of this paper). In the Ref. [8] 
the Eqs. (4) arc solved. This solution is known since Filer's time. The essence of the problem 
is in the description of the dynamic of the process, i.e., in approximation (12). 

In the abstract [8] we read:" ...within a simple model which respects unitarity". Unitarity 
within a model with " anti-Hermitian" operator Un — — iF/2? 

Also Cal writes: "Dover, Gal, Richard [9] pinpointed errors in Ref. [4]." In fact, the situation 
is inverse. The simple calculation given in [9] confirms the FTA and was already done in our 
first paper [4] (zero angle rescattering diagrams of n) as well as in Ref. [5] (see Pgs. R1884, 
R1885). Since the authors [8,9] "overlooked" all of this we cite here only one paragraph from 
Ref. [5]: "The authors [9] substitute H = -iF/2 in Eq.(4) of Ref. [5] and obtain Wpot- On the 
basis of this and only this they refute the result of Ref. [4]. In other words they refute our result 
because it differs from the potential model one." 

Unusual logic. The procedure mentioned above is a verification of FTA by the example of 
exactly solveble (but incorrect) potential model. We abandon the approximation (12), i.e., the 
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potential description of n-medium interaction in principal, what was clearly pointed in [4,5]. 

8 On an anhancement of oscillations 

In the region Ft ^ 1 

Wann{t)/Wj,ot{t) ~ -ImVt ~ H > 1 

what means an anhancement of oscillations as compared with standard result. Certainly, this 
is a record value. For other problems an enhancement factor will provoke a lesser skepticism. 
In fact the suppression factor Rpot ~ 10~^^ is surprising. 

When Ft < 1, annihilation products as well as n in the final state can be detected. In this 
case condition (24) is not fulfilled. Nevertheless, the qualitative picture remains as before: in 
the standard approach ImUn leads to suppression of oscillations, in our one this is not the case. 
The presence of open channels of n-medium interaction (which are described by ImUn in the 
potential model) do not lead to suppression of oscillations. 

For large part of actual problems the both approaches give an identical result. For two-step 
processes of the type (1), when 6-particle absorption is essential, the principal disagreement 
takes place. We do not want to make the categorical conclusions. Nevertheless, two points are 
obvious: (1) The potential approach should be considered only as starting point. (2) The only 
assumption used in our calculation is expression (24). However, this condition is basic to the 
all (S-matrix calculations. 
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Figure caption 

Fig.l. (a) nn transition in a medium followed by annihilation; (b) The same as in Fig. (a), but 
with escaping of particle in the nn transition vertex. 
Fig. 2. Free-space process n ^ n ^ pe^v. 
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